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I Abstract. The A-scaling method has been applied to the total multiplicity 

distribution of the relativistic ion collisions of p+p, C+C and Pb+Pb which were 
simulated by a Monte Carlo package, LUCIAE 3.0. It is found that the A-scaling 
parameter decreases with the increasing of the system size. Moreover, the heat 
. capacities of different mesons and baryons have been extracted from the event-by- 

event temperature fluctuation in the region of low transverse mass and they show the 
dropping trend with the increasing of impact parameter. 
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^ I Studying fluctuation and correlation is an important issue to investigate matter 

^ properties formed in relativistic nucleus-nucleus collisions. Some observables have 

been suggested to quantify the fluctuation, such as the balance function, net charge 
^ ' fluctuation, multiplicity fluctuation, transverse momentum fluctuation, particle ratio 
' fluctuation and $pt etc pLj. In this article we will discuss some novel methods. We will 
apply a universal fluctuation method (A-scaling) which was recently proposed to explore 
the order- disorder phase transition in the low-intermediate energy heavy ion collisions 
[3 El to discuss multiplicity fluctuation in relativistic ion collisions. In addition, the heat 
capacity of hadrons are extracted from the event-by-event temperature fluctuation. 

A-scaling is observed when two or more probability distributions P[m] of the 
stochastic observable m collapse onto a single scaling curve $(z) if a new scaling 
observable is defined by 2; = -^^^p-^. This curve is 0: 

(m) V[m] = ^ $[^^^] (1) 

where A is a scaling parameter, m* is the most probable value of m, and (m) 
is the mean of m. If we assume that P[m] is a Gaussian distribution, we have 
P[m] = g.^/2^ 6^P[~|(^^^v^)^]i where /i = (m) = m*, a is the width of the Gaussian 
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Figure 1. The determination of the A-scaling parameter using the power-law fits to 
the relationship mjiMS vs (m). From the left to right, it corresponds to p + p with 
A = 1.35, C + C with A = 1.00 and Pb+Pb with A = 0.80. The circles are the 
calculated values and the lines represent the power-law fits. 

distribution, both depending on incident energy. If this Gaussian distribution P[m] 
obeys A-scahng law, we have 

a oc /i^. (2) 

In this work, p+p, C+C and Pb-hPb in SPS energies (20-200 AGeV) have been 
investigated with the help of LUCIAE3.0 model The head-on (b=Ofm) collisions 
are simulated in this work. The LUCIAE is an extension Monte Carlo model of the 
FRITIOF jHj, in which a nucleus - nucleus collision is described as the independent sum 
of nucleon-nucleon collisions. For more details of LUCIAE, please see Ref. jl]. 

In order to better investigate the A-scaling, we can use the power-law fit to the 
relationship between the RMS width of the multiplicity distribution {mjiMs) and the 
mean value of the multiplicity distribution ((m)) according to Eq.(l). If there exists a 
good fit: rriRMS = co{m)^ where cq is a fit parameter, indicating that a good A-scaling 
law is satisfied with a scaling parameter A. In this way, we found that the best fit value 
of A is 1.35 for p-|-p, 1.00 for C-l-C and 0.80 for Pb -|- Pb, respectively, which is shown 
in Figure 1. In the figure, each point represents an energy point from 20 AGeV to 200 
AGeV with the interval of 20 AGeV. When A = 1, the RMS width (fluctuation) is 
proportional to its mean value. This is the case of C-l-C system. For light system p+p, 
the rising of the RMS width is faster than its mean value with the increasing of beam 
energy. In contrary, for heavier system, the growing of the RMS width is slower than 
the mean value. In all three cases. A- value is larger than 1/2, which is the case of the 
independent particle emission (Poisson distribution). The decreasing behaviour of A- 
value with the system size reflects how the multiplicity fluctuation grows with its mean 
value. Its physics origin may stem from the stronger rescattering effect in larger system 
which induce some correlations in multiparticle production process and in consequence 
the fluctuation becomes weaken. 

Using these values of A, the multiplicity distributions can be compressed onto a 
unique curve for a given collision system from 20 AGeV to 200 AGeV. Figure 2 shows 
the A-scaling curves for p+p, C-l-C and Pb-|-Pb systems with A value of 1.35, 1.00 and 
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0.80, respectively. 

In relativistic collisions, the transverse mass (m^") distribution of the particles in 
lower rriT region can be approximately described as — = Ae~~ where A is a 
normalized coefficient which is related to the volume/multiplicity term and T is an 
apparent temperature. Thanks to the large enough particle multiplicity, it is possible 
to extract A and T for particles on event-by-event basis. In this way, the event-by- 
event temperature distribution can be constructed. Usually, such kind of temperature 
distribution (P(T)) can be described by [HI EI 

AT 

P(T) ~ exp[-C,{—n (3) 

where AT is the deviation (T — (T)) of temperature from the mean value ((T)) and 
can be explained as the heat capacity of a certain hadron. Since the heat capacity is an 
extensive observable, we define a normalized heat capacity C^/N, i.e. the heat capacity 
per hadron multiplicity (A^) which corresponds to the energy that is required to go up 
one unit temperature for producing one particle. 

The heat capacity of different particles (in the range of rriT — rriQ < 0.5 GeV) has 
been extracted by Eq.(3). The left panel of Figure 3 shows an example of event-by-event 
temperature distribution of vr"*" for the head-on Pb+Pb collisions at 160 AGeV. From 
this kind of distribution, we apply Eq.(3) to extract the heat capacity C^. The curve of 
the figure depicts this fit. In this way, we get the heat capacity for various hadrons in 
different energies. The middle panel of Figure 3 gives the dependence of heat capacity 
of various particles (tt"'", k~^, p and A) on incident energy in head-on Pb-fPb collision. 
From this figure, the normalized heat capacity of various particles decreases with the 
increasing of incident energy and tends to a saturation, indicating that less energy is 
needed to rise the same temperature in higher energies. In the viewpoint of fluctuation, 
this means that the increasing ttlt fluctuation at higher energies. The normalized heat 
capacity of the particles increases with mass of particle, which reflects that more energy 
is needed to rise the same temperature for heavier particles. In addition, the impact 
parameter dependence of the heat capacity is also investigated for vr'*'. The right panel 
of Figure 3 depicts the dependence oi C^/N for tt"*" on the impact parameter for Pb + 
Pb. The decreasing trend of C^/N indicates that the temperature fluctuation becomes 
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Figure 3. (a) The event-by-event temperature distribution of 7r+ in the head-on 
Pb-|-Pb coUision at 160AGeV. The curve is a fit with Eq.(3); (b) The normahzed heat 
capacity of different mesons and baryons as a function of incident energy for the head- 
on Pb + Pb coUisions. The hues represent the second polynomial fits to guide the 
eyes; (c) The dependence of the normalized heat capacity of 7r+ on impact parameter 
for 160 AGeV Pb + Pb collisions. The error bars are statistical. 



larger in peripheral collisions, which is similar to the Pt fluctuation or balance function 
as a function of centrality in STAR data [S]. 

In summary, we have demonstrated that there exists the A-scaling for the total 
multiplicity of charged particles for the simulated head-on collisions of p + p, C + 
C and Pb + Pb from Euh = 20 to 200 AGeV in LUCIAE model. That means 
that the multiplicity distributions obey a certain kind of universal laws, regardless 
of beam energies and collision systems. It is found that the scaling value of the A 
decreases with the system size, which reflects that the growth of the fluctuation with 
the multiplicity is faster in light system. The event-by-event thermal fluctuation in lower 
rriT region is constructed and the normalized heat capacities of the different mesons and 
baryons are extracted. It is found that the heat capacity per hadron decreases with 
the increasing incident energy and impact parameter, while it increases with mass of 
particle. Considering that the LUCIAE model is a string-hadronic model, the partonic 
effect which becomes more important in RHIC energies is absent in the model. This 
effect should be investigated in near future. The work along this line is in progress. 
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